Evaluation of specific RBE in different cells of hippocampus under high-dose proton irradiation in rats.

The study aimed to determine the specific relative biological effectiveness (RBE) of various cells in the hippocampus following proton irradiation. Sixty Sprague-Dawley rats were randomly allocated to 5 groups receiving 20 or 30 Gy of proton or photon irradiation. Pathomorphological neuronal damage in the hippocampus was assessed using Hematoxylin-eosin (HE) staining. The expression level of NeuN, Nestin, Caspase-3, Olig2, CD68 and CD45 were determined by immunohistochemistry (IHC). The RBE range established by comparing the effects of proton and photon irradiation at equivalent biological outcomes. Proton20Gy induced more severe damage to neurons than photon20Gy, but showed no difference compared to photon30Gy. The RBE of neuron was determined to be 1.65. Similarly, both proton20Gy and proton30Gy resulted in more inhibition of oligodendrocytes and activation of microglia in the hippocampal regions than photon20Gy and photon30Gy. However, the expression of Olig2 was higher and CD68 was lower in the proton20Gy group than in the photon30Gy group. The RBE of oligodendrocyte and microglia was estimated to be between 1.1 to 1.65. For neural stem cells (NSCs) and immune cells, there were no significant difference in the expression of Nestin and CD45 between proton and photon irradiation (both 20 and 30 Gy). Therefore, the RBE for NSCs and immune cell was determined to be 1.1. These findings highlight the varying RBE values of different cells in the hippocampus in vivo. Moreover, the actual RBE of the hippocampus may be higher than 1.1, suggesting that using as RBE value of 1.1 in clinical practice may underestimate the toxicities induced by proton radiation.

Benchmarking proton RBE models.

Objective.To biologically optimise proton therapy, models which can accurately predict variations in proton relative biological effectiveness (RBE) are essential. Current phenomenological models show large disagreements in RBE predictions, due to different model assumptions and differences in the data to which they were fit. In this work, thirteen RBE models were benchmarked against a comprehensive proton RBE dataset to evaluate predictions when all models are fit using the same data and fitting techniques, and to assess the statistical robustness of the models.Approach.Model performance was initially evaluated by fitting to the full dataset, and then a cross-validation approach was applied to assess model generalisability and robustness. The impact of weighting the fit and the choice of biological endpoint (either single or multiple survival levels) was also evaluated.Main results.Fitting the models to a common dataset reduced differences between their predictions, however significant disagreements remained due to different underlying assumptions. All models performed poorly under cross-validation in the weighted fits, suggesting that some uncertainties on the experimental data were significantly underestimated, resulting in over-fitting and poor performance on unseen data. The simplest model, which depends linearly on the LET but has no tissue or dose dependence, performed best for a single survival level. However, when fitting to multiple survival levels simultaneously, more complex models with tissue dependence performed better. All models had significant residual uncertainty in their predictions compared to experimental data.Significance.This analysis highlights that poor quality of error estimation on the dose response parameters introduces substantial uncertainty in model fitting. The significant residual error present in all approaches illustrates the challenges inherent in fitting to large, heterogeneous datasets and the importance of robust statistical validation of RBE models.

Dose and dose rate dependence of the tissue sparing effect at ultra-high dose rate studied for proton and electron beams using the zebrafish embryo model.

PURPOSE: A better characterization of the dependence of the tissue sparing effect at ultra-high dose rate (UHDR) on physical beam parameters (dose, dose rate, radiation quality) would be helpful towards a mechanistic understanding of the FLASH effect and for its broader clinical translation. To address this, a comprehensive study on the normal tissue sparing at UHDR using the zebrafish embryo (ZFE) model was conducted. METHODS: One-day-old ZFE were irradiated over a wide dose range (15-95 Gy) in three different beams (proton entrance channel, proton spread out Bragg peak and 30 MeV electrons) at UHDR and reference dose rate. After irradiation the ZFE were incubated for 4 days and then analyzed for four different biological endpoints (pericardial edema, curved spine, embryo length and eye diameter). RESULTS: Dose-effect curves were obtained and a sparing effect at UHDR was observed for all three beams. It was demonstrated that proton relative biological effectiveness and UHDR sparing are both relevant to predict the resulting dose response. Dose dependent FLASH modifying factors (FMF) for ZFE were found to be compatible with rodent data from the literature. It was found that the UHDR sparing effect saturates at doses above âˆ¼ 50 Gy with an FMF of âˆ¼ 0.7-0.8. A strong dose rate dependence of the tissue sparing effect in ZFE was observed. The magnitude of the maximum sparing effect was comparable for all studied biological endpoints. CONCLUSION: The ZFE model was shown to be a suitable pre-clinical high-throughput model for radiobiological studies on FLASH radiotherapy, providing results comparable to rodent models. This underlines the relevance of ZFE studies for FLASH radiotherapy research.

Beam quality correction factors for ionization chambers in a 0.35 T magnetic resonance (MR)-linac - A Monte Carlo study.

PURPOSE: The purpose of this study was to directly calculate [Formula: see text] correction factors for four cylindrical ICs for a 0.35 T MR-linac using the Monte Carlo (MC) method. METHODS: A previously-validated TOPAS/GEANT4 MC head model of the 0.35 T MR-linac was employed. The MR-compatible Exradin A12, A1SL, A26, and A28 cylindrical ICs were modeled considering the dead volume in the air cavity. The [Formula: see text] correction factor was determined for initial electron energies of 5-7 MeV. The correction factor was calculated for all four angular orientations in the lateral plane. The impact of the 0.35 T magnetic field on the IC response was also investigated. RESULTS: The maximum beam quality dependence in the [Formula: see text] exhibited by the A12, A1SL, A26, and A28 ICs was 1.10 %, 2.17 %, 0.81 %, and 1.75 %, respectively, considering all angular orientations. The magnetic field dependence was < 1 % and the maximum [Formula: see text] correction was < 2 % when the detector was aligned along the direction of the magnetic field at 0° and 180° angles. The A12 IC over-responded up to 5.40 % for the orthogonal orientation. An asymmetry in the response of up to 8.30 % was noted for the A28 IC aligned at 90° and 270° angles. CONCLUSIONS: A parallel orientation for the IC, with respect to the magnetic field, is recommended for reference dosimetry in MRgRT. Both over and under-response in the IC signal was noted for the orthogonal orientations, which is highly dependent on the cavity diameter, cavity length, and the dead volume.

Equivalent uniform RBE-weighted dose in eye plaque brachytherapy.

BACKGROUND: Brachytherapy for ocular melanoma is based on the application of eye plaques with different spatial dose nonuniformity, time-dependent dose rates and relative biological effectiveness (RBE). PURPOSE: We propose a parameter called the equivalent uniform RBE-weighted dose (EUDRBE) that can be used for quantitative characterization of integrated cell survival in radiotherapy modalities with the variable RBE, dose nonuniformity and dose rate. The EUDRBE is applied to brachytherapy with 125I eye plaques designed by the Collaborative Ocular Melanoma Study (COMS). METHODS: The EUDRBE is defined as the uniform dose distribution with RBE = 1 that causes equal cell survival for a given nonuniform dose distribution with the variable RBE > 1. The EUDRBE can be used for comparison of cell survival for nonuniform dose distributions with different RBE, because they are compared to the reference dose with RBE = 1. The EUDRBE is applied to brachytherapy with 125I COMS eye plaques that are characterized by a steep dose gradient in tumor base-apex direction, protracted irradiation during time intervals of 3-8 days, and variable dose-rate dependent RBE with a maximum of about 1.4. The simulations are based on dose of 85 Gy prescribed to the farthest intraocular extent of the tumor (tumor apex). To compute the EUDRBE in eye plaque brachytherapy and correct for protracted irradiation, the distributions of physical dose have been converted to non-uniform distributions of biologically effective dose (BED) to include the biological effects of sublethal cellular repair, Our radiobiological analysis considers the combined effects of different time-dependent dose rates, spatial dose non-uniformity, dose fractionation and different RBE and can be used to derive optimized dose regimens brachytherapy. RESULTS: Our simulations show that the EUDRBE increases with the prescription depths and the maximum increase may achieve 6% for the tumor height of 12 mm. This effect stems from a steep dose gradient within the tumor that increases with the prescription depth. The simulations also show that the EUDRBE increase may achieve 12% with increasing the dose rate when implant duration decreases. The combined effect of dose nonuniformity and dose rate may change the EUDRBE up to 18% for the same dose prescription of 85 Gy to tumor apex. The absolute dose range of 48-61 Gy (RBE) for the EUDRBE computed using 4 or 5 fractions is comparable to the dose prescriptions used in stereotactic body radiation therapy (SBRT) with megavoltage X-rays (RBE = 1) for different cancers. The tumor control probabilities in SBRT and eye plaque brachytherapy are very similar at the level of 80% or higher that support the hypothesis that the selected approximations for the EUDRBE are valid. CONCLUSIONS: The computed range of the EUDRBE in 125I COMS eye plaque brachytherapy suggests that the selected models and hypotheses are acceptable. The EUDRBE can be useful for analysis of treatment outcomes and comparison of different dose regimens in eye plaque brachytherapy.

Primary study of the relative and compound biological effectiveness model for boron neutron capture therapy based on nanodosimetry.

BACKGROUND: The current radiobiological model employed for boron neutron capture therapy (BNCT) treatment planning, which relies on microdosimetry, fails to provide an accurate representation the biological effects of BNCT. The precision in calculating the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) plays a pivotal role in determining the therapeutic efficacy of BNCT. Therefore, this study focuses on how to improve the accuracy of the biological effects of BNCT. PURPOSE: The purpose of this study is to propose new radiation biology models based on nanodosimetry to accurately assess RBE and CBE for BNCT. METHODS: Nanodosimetry, rooted in ionization cluster size distributions (ICSD), introduces a novel approach to characterize radiation quality by effectively delineating RBE through the ion track structure at the nanoscale. In the context of prior research, this study presents a computational model for the nanoscale assessment of RBE and CBE. We establish a simplified model of DNA chromatin fiber using the Monte Carlo code TOPAS-nBio to evaluate the applicability of ICSD to BNCT and compute nanodosimetric parameters. RESULTS: Our investigation reveals that both homogeneous and heterogeneous nanodosimetric parameters, as well as the corresponding biological model coefficients α and ß, along with RBE values, exhibit variations in response to varying intracellular 10B concentrations. Notably, the nanodosimetric parameter M 1 C 2 $M_1^{{{\mathrm{C}}}_2}$ effectively captures the fluctuations in model coefficients α and RBE. CONCLUSION: Our model facilitates a nanoscale analysis of BNCT, enabling predictions of nanodosimetric quantities for secondary ions as well as RBE, CBE, and other essential biological metrics related to the distribution of boron. This contribution significantly enhances the precision of RBE calculations and holds substantial promise for future applications in treatment planning.

Hypoxia-informed RBE-weighted beam orientation optimization for intensity modulated proton therapy.

BACKGROUND: Variable relative biological effectiveness (RBE) models in treatment planning have been proposed to optimize the therapeutic ratio of proton therapy. It has been reported that proton RBE decreases with increasing tumor oxygen level, offering an opportunity to address hypoxia-related radioresistance with RBE-weighted optimization. PURPOSE: Here, we obtain a voxel-level estimation of partial oxygen pressure to weigh RBE values in a single biologically informed beam orientation optimization (BOO) algorithm. METHODS: Three glioblastoma patients with [18 F]-fluoromisonidazole (FMISO)-PET/CT images were selected from the institutional database. Oxygen values were derived from tracer uptake using a nonlinear least squares curve fitting. McNamara RBE, calculated from proton dose, was then weighed using oxygen enhancement ratios (OER) for each voxel and incorporated into the dose fidelity term of the BOO algorithm. The nonlinear optimization problem was solved using a split-Bregman approach, with FISTA as the solver. The proposed hypoxia informed RBE-weighted method (HypRBE) was compared to dose fidelity terms using the constant RBE of 1.1 (cRBE) and the normoxic McNamara RBE model (RegRBE). Tumor homogeneity index (HI), maximum biological dose (Dmax), and D95%, as well as OAR therapeutic index (TI = gEUDCTV /gEUDOAR ) were evaluated along with worst-case statistics after normalization to normal tissue isotoxicity. RESULTS: Compared to [cRBE, RegRBE], HypRBE increased tumor HI, Dmax, and D95% across all plans by on average [31.3%, 31.8%], [48.6%, 27.1%], and [50.4%, 23.8%], respectively. In the worst-case scenario, the parameters increase on average by [12.5%, 14.7%], [7.3%,-8.9%], and [22.3%, 2.1%]. Despite increased OAR Dmean and Dmax by [8.0%, 3.0%] and [13.1%, -0.1%], HypRBE increased average TI by [22.0%, 21.1%]. Worst-case OAR Dmean, Dmax, and TI worsened by [17.9%, 4.3%], [24.5%, -1.2%], and [9.6%, 10.5%], but in the best cases, HypRBE escalates tumor coverage significantly without compromising OAR dose, increasing the therapeutic ratio. CONCLUSIONS: We have developed an optimization algorithm whose dose fidelity term accounts for hypoxia-informed RBE values. We have shown that HypRBE selects bE:\Alok\aaeams better suited to deliver high physical dose to low RBE, hypoxic tumor regions while sparing the radiosensitive normal tissue.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool.

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Effects of cellular radioresponse on therapeutic helium-, carbon-, oxygen-, and neon-ion beams: a simulation study.

Objective. Helium, oxygen, and neon ions in addition to carbon ions will be used for hypofractionated multi-ion therapy to maximize the therapeutic effectiveness of charged-particle therapy. To use new ions in cancer treatments based on the dose-fractionation protocols established in carbon-ion therapy, this study examined the cell-line-specific radioresponse to therapeutic helium-, oxygen-, and neon-ion beams within wide dose ranges.Approach. Response of cells to ions was described by the stochastic microdosimetric kinetic model. First, simulations were made for the irradiation of one-field spread-out Bragg peak beams in water with helium, carbon, oxygen, and neon ions to achieve uniform survival fractions at 37%, 10%, and 1% for human salivary gland tumor (HSG) cells, the reference cell line for the Japanese relative biological effectiveness weighted dose system, within the target region defined at depths from 90 to 150 mm. The HSG cells were then replaced by other cell lines with different radioresponses to evaluate differences in the biological dose distributions of each ion beam with respect to those of carbon-ion beams.Main results. For oxygen- and neon-ion beams, the biological dose distributions within the target region were almost equivalent to those of carbon-ion beams, differing by less than 5% in most cases. In contrast, for helium-ion beams, the biological dose distributions within the target region were largely different from those of carbon-ion beams, more than 10% in several cases.Significance.From the standpoint of tumor control evaluated by the clonogenic cell survival, this study suggests that the dose-fractionation protocols established in carbon-ion therapy could be reasonably applied to oxygen- and neon-ion beams while some modifications in dose prescription would be needed when the protocols are applied to helium-ion beams. This study bridges the gap between carbon-ion therapy and hypofractionated multi-ion therapy.

Cell Survival Computation via the Generalized Stochastic Microdosimetric Model (GSM2); Part II: Numerical Results.

In the present paper we numerically investigate, using Monte Carlo simulation, the theoretical results predicted by the Generalized Stochastic Microdosimetric Model (GSM2), as shown in the published companion paper. Taking advantage of the particle irradiation data ensemble (PIDE) dataset, we calculated GSM2 biological parameters of human salivary gland (HSG) and V79 cell lines. Further, exploiting the TOPAS-microdosimetric extension, we simulated the microdosimetric spectra of different radiation fields of therapeutic interest generated by four different ions (protons, helium-4, carbon-12 and oxygen-16) each at three different residual ranges. We investigated the properties of the initial damage distributions as well as the cell survival curve predicted by GSM2, focusing especially on the non-Poissonian effects naturally included in the model. GSM2 successfully computed cell survival curves, accurately describing experimental behavior even under challenging LET and dose conditions.

Enhanced RBE of Particle Radiation Depends on Beam Size in the Micrometer Range.

High-linear energy transfer (LET) radiation, such as heavy ions is associated with a higher relative biological effectiveness (RBE) than low-LET radiation, such as photons. Irradiation with low- and high-LET particles differ in the interaction with the cellular matter and therefore in the spatial dose distribution. When a single high-LET particle interacts with matter, it results in doses of up to thousands of gray (Gy) locally concentrated around the ion trajectory, whereas the mean dose averaged over the target, such as a cell nucleus is only in the range of a Gy. DNA damage therefore accumulates in this small volume. In contrast, up to hundreds of low-LET particle hits are required to achieve the same mean dose, resulting in a quasi-homogeneous damage distribution throughout the cell nucleus. In this study, we investigated the dependence of RBE from different spatial dose depositions using different focused beam spot sizes of proton radiation with respect to the induction of chromosome aberrations and clonogenic cell survival. Human-hamster hybrid (AL) as well as Chinese hamster ovary cells (CHO-K1) were irradiated with focused low LET protons of 20 MeV (LET = 2.6 keV/µm) beam energy with a mean dose of 1.7 Gy in a quadratic matrix pattern with point spacing of 5.4 × 5.4 µm2 and 117 protons per matrix point at the ion microbeam SNAKE using different beam spot sizes between 0.8 µm and 2.8 µm (full width at half maximum). The dose-response curves of X-ray reference radiation were used to determine the RBE after a 1.7 Gy dose of radiation. The RBE for the induction of dicentric chromosomes and cell inactivation was increased after irradiation with the smallest beam spot diameter (0.8 µm for chromosome aberration experiments and 1.0 µm for cell survival experiments) compared to homogeneous proton radiation but was still below the RBE of a corresponding high LET single ion hit. By increasing the spot size to 1.6-1.8 µm, the RBE decreased but was still higher than for homogeneously distributed protons. By further increasing the spot size to 2.7-2.8 µm, the RBE was no longer different from the homogeneous radiation. Our experiments demonstrate that varying spot size of low-LET radiation gradually modifies the RBE. This underlines that a substantial fraction of enhanced RBE originates from inhomogeneous energy concentrations on the µm scale (mean intertrack distances of low-LET particles below 0.1 µm) and quantifies the link between such energy concentration and RBE. The missing fraction of RBE enhancement when comparing with high-LET ions is attributed to the high inner track energy deposition on the nanometer scale. The results are compared with model results of PARTRAC and LEM for chromosomal aberration and cell survival, respectively, which suggest mechanistic interpretations of the observed radiation effects.

Sensitivity of a mini-TEPC to radiation quality variations in clinical proton beams.

PURPOSE: This work aims at studying the sensitivity of a miniaturized Tissue-Equivalent Proportional Counter to variations of beam quality in clinical radiation fields, to further investigate its performances as radiation quality monitor. METHODS: Measurements were taken at the CATANA facility (INFN-LNS, Catania, Italy), in a monoenergetic and an energy-modulated proton beam with the same initial energy of 62 MeV. PMMA layers were placed in front of the detector to measure at different depths along the depth-dose profile. The frequency- and dose-mean lineal energy were compared to the track- and dose-averaged LET calculated by Monte Carlo simulations. A microdosimetric evaluation of the Relative Biological Effectiveness (RBE) was performed and compared with cell survival experiments. RESULTS: Microdosimetric distributions measured at identical depths in the two beams show spectral differences reflecting their different radiation quality. Discrepancies are most evident at depths corresponding to the Spread-Out Bragg Peak, while spectra at the entrance and in the dose fall-off regions are similar. This can be explained by the different energy components that compose the pristine and spread-out peaks at each depth. The trend of microdosimetric mean values matches that of calculated LET averages along the entire penetration depth, and the microdosimetric estimation of RBE is consistent with radiobiological data not only at 2 Gy but also at lower dose levels, such as those absorbed by healthy tissues. CONCLUSIONS: The mini-TEPC is sensitive to differences in radiation quality resulting from different modulations of the proton beam, confirming its potential for beam quality monitoring in proton therapy.

Proton relative biological effectiveness for the induction of DNA double strand breaks based on Geant4.

Uncertainties in the relative biological effectiveness (RBE) of proton remains a major barrier to the biological optimization of proton therapy. A large amount of experimental data suggest that proton RBE is variable. As an evolving Monte Carlo code toolkit, Geant4-DNA is able to simulate the initial DNA damage caused by particle beams through physical and chemical interactions at the nanometer scale over a short period of time. This contributes to evaluating the radiobiological effects induced by ionizing radiation. Based on the Geant4-DNA toolkit, this study constructed a DNA geometric model containing 6.32Gbp, simulated the relationship between radiochemical yields (G-values) and their corresponding chemical constructors, and calculated a detailed calculation of the sources of damage and the complexity of damage in DNA strand breaks. The damage model constructed in this study can simulate the relative biological effectiveness (RBE) in the proton Bragg peak region. The results indicate that: (1) When the electron energy is below 400 keV, the yield of OH·account for 18.1% to 25.3% of the total water radiolysis yields. (2) Under the influence of histone clearance function, the yield of indirect damage account for over 72.93% of the yield of DNA strand breaks (SBs). When linear energy transfer (LET) increased from 29.79 (keV/µm) to 64.29 (keV/µm), the yield of double strand breaks (DSB) increased from 17.27% to 32.65%. (3) By investigating the effect of proton Bragg peak depth on the yield of direct DSB (DSBdirect) and total DSB (DSBtotal), theRBEDSBtotandRBEDSBdirlevels of cells show that the RBE value of protons reaches 2.2 in the Bragg peak region.

Boron-Containing MOF Nanoparticles with Stable Metabolism in U87-MG Cells Combining Microdosimetry To Evaluate Relative Biological Effectiveness of Boron Neutron Capture Therapy.

Accurate prediction of the relative biological effectiveness (RBE) of boron neutron capture therapy (BNCT) is challenging. The therapy is different from other radiotherapy; the dynamic distribution of boron-containing compounds in tumor cells affects the therapeutic outcome considerably and hampers accurate measurement of the neutron-absorbed dose. Herein, we used boron-containing metal-organic framework nanoparticles (BMOFs) with high boron content to target U87-MG cells and maintain the concentration of the 10B isotope in cells. The content of boron in the cells could maintain 90% (60 ppm) within 20 min compared with that at the beginning; therefore, the accurate RBE of BNCT can be acquired. The effects of BNCT upon cells after neutron irradiation were observed, and the neutron-absorbed dose was obtained by Monte Carlo simulations. The RBE of BMOFs was 6.78, which was 4.1-fold higher than that of a small-molecule boron-containing agent (boric acid). The energy spectrum of various particles was analyzed by Monte Carlo simulations, and the RBE was verified theoretically. Our results suggested that the use of nanoparticle-based boron carriers in BNCT may have many advantages and that maintaining a stable boron distribution within cells may significantly improve the efficiency of BNCT.

Generalized methods for predicting biological response to mixed radiation types and calculating equieffective doses (EQDX).

BACKGROUND: Predicting biological responses to mixed radiation types is of considerable importance when combining radiation therapies that use multiple radiation types and delivery regimens. These may include the use of both low- and high-linear energy transfer (LET) radiations. A number of theoretical models have been developed to address this issue. However, model predictions do not consistently match published experimental data for mixed radiation exposures. Furthermore, the models are often computationally intensive. Accordingly, there is a need for efficient analytical models that can predict responses to mixtures of low- and high-LET radiations. Additionally, a general formalism to calculate equieffective dose (EQDX) for mixed radiations is needed. PURPOSE: To develop a computationally efficient analytical model that can predict responses to complex mixtures of low- and high-LET radiations as a function of either absorbed dose or EQDX. METHODS: The Zaider-Rossi model (ZRM) was modified by replacing the geometric mean of the quadratic coefficients in the interaction term with the arithmetic mean. This modified ZRM model (mZRM) was then further generalized to any number of radiation types and its validity was tested against published experimental observations. Comparisons between the predictions of the ZRM and mZRM, and other models, were made using two and three radiation types. In addition, a generalized formalism for calculating EQDX for mixed radiations was developed within the context of mZRM and validated with published experimental results. RESULTS: The predictions of biological responses to mixed-LET radiations calculated with the mZRM are in better agreement with experimental observations than ZRM, especially when high- and low-LET radiations are mixed. In these situations, the ZRM overestimated the surviving fraction. Furthermore, the EQDX calculated with mZRM are in better agreement with experimental observations. CONCLUSION: The mZRM is a computationally efficient model that can be used to predict biological response to mixed radiations that have low- and high-LET characteristics. Importantly, interaction terms are retained in the calculation of EQDX for mixed radiation exposures within the mZRM framework. The mZRM has application in a wide range of radiation therapies, including radiopharmaceutical therapy.

Computational approaches in the estimation of radiobiological damage for human-malignant cells irradiated with clinical proton and carbon beams.

PURPOSE: The use of Monte Carlo (MC) simulations capable of reproducing radiobiological effects of ionising radiation on human cell lines is of great importance, especially for cases involving protons and heavier ion beams. In the latter, huge uncertainties can arise mainly related to the effects of the secondary particles produced in the beam-tissue interaction. This paper reports on a detailed MC study performed using Geant4-based approach on three cancer cell lines, the HTB-177, CRL-5876 and MCF-7, that were previously irradiated with therapeutic proton and carbon ion beams. METHODS: A Geant4-based approach used jointly with analytical calculations has been developed to provide a more realistic estimation of the radiobiological damage produced by proton and carbon beams in tissues, reproducing available data obtained from in vitro cell irradiations. The MC "Hadrontherapy" Geant4 application and the Local Effect Model: LEM I, LEM II and LEM III coupled with the different numerical approaches: RapidRusso (RR) and RapidScholz (RS) were used in the study. RESULTS: Experimental survival curves are compared with those evaluated using the highlighted Geant4 MC-based approach via chi-square statistical analysis, for the combinations of radiobiological models and numerical approaches, as outlined above. CONCLUSION: This study has presented a comparison of the survival data from MC simulations to experimental survival data for three cancer cell lines. An overall best level of agreement was obtained for the HTB-177 cells.

Monte Carlo calculated beam quality correction factors for high energy electron beams.

PURPOSE: As the dosimetry protocol TRS 398 is being revised and the ICRU report 90 provides new recommendations for density correction as well as the mean ionization energies of water and graphite, updated beam quality correction factors kQ are calculated for reference dosimetry in electron beams and for independent validation of previously determined values. METHODS: Monte Carlo simulations have been performed using EGSnrc to calculate the absorbed dose to water and the dose to the active volumes of ionization chambers SNC600c, SNC125c and SNC350p (all Sun Nuclear, A Mirion Medical Company, Melbourne, FL). Realistic clinical electron beam spectra were used to cover the entire energy range of therapeutic electron accelerators. The Monte Carlo simulations were validated by measurements on a clinical linear accelerator. With regards to the cylindrical chambers, the simulations were performed according to the setup recommendations of TRS 398 and AAPM TG 51, i.e. with and without consideration of a reference point shift by rcav/2. RESULTS: kQ values as a function of the respective beam quality specifier R50 were fitted by recommended equations for electron beam dosimetry in the range of 5 MeV to 18 MeV. The fitting curves to the calculated values showed a root mean square deviation between 0.0016 and 0.0024. CONCLUSION: Electron beam quality correction factors kQ were calculated by Monte Carlo simulations for the cylindrical ionization chambers SNC600c and SNC125c as well as the plane parallel ionization chamber SNC350p to provide updated data for the TRS 398 and TG 51 dosimetry protocols.

Incorporating variable RBE in IMPT optimization for ependymoma.

PURPOSE: To study the dosimetric impact of incorporating variable relative biological effectiveness (RBE) of protons in optimizing intensity-modulated proton therapy (IMPT) treatment plans and to compare it with conventional constant RBE optimization and linear energy transfer (LET)-based optimization. METHODS: This study included 10 pediatric ependymoma patients with challenging anatomical features for treatment planning. Four plans were generated for each patient according to different optimization strategies: (1) constant RBE optimization (ConstRBEopt) considering standard-of-care dose requirements; (2) LET optimization (LETopt) using a composite cost function simultaneously optimizing dose-averaged LET (LETd ) and dose; (3) variable RBE optimization (VarRBEopt) using a recent phenomenological RBE model developed by McNamara et al.; and (4) hybrid RBE optimization (hRBEopt) assuming constant RBE for the target and variable RBE for organs at risk. By normalizing each plan to obtain the same target coverage in either constant or variable RBE, we compared dose, LETd , LET-weighted dose, and equivalent uniform dose between the different optimization approaches. RESULTS: We found that the LETopt plans consistently achieved increased LET in tumor targets and similar or decreased LET in critical organs compared to other plans. On average, the VarRBEopt plans achieved lower mean and maximum doses with both constant and variable RBE in the brainstem and spinal cord for all 10 patients. To compensate for the underdosing of targets with 1.1 RBE for the VarRBEopt plans, the hRBEopt plans achieved higher physical dose in targets and reduced mean and especially maximum variable RBE doses compared to the ConstRBEopt and LETopt plans. CONCLUSION: We demonstrated the feasibility of directly incorporating variable RBE models in IMPT optimization. A hybrid RBE optimization strategy showed potential for clinical implementation by maintaining all current dose limits and reducing the incidence of high RBE in critical normal tissues in ependymoma patients.

Impact of gadolinium concentration and cell oxygen levels on radiobiological characteristics of gadolinium neutron capture therapy technique in brain tumor treatment.

Neutron capture therapy (NCT) with various concentrations of gadolinium (157Gd) is one of the treatment modalities for glioblastoma (GBM) tumors. Current study aims to evaluate how variations of 157Gd concentration and cell oxygen levels can affect the relative biological effectiveness (RBE) of gadolinium neutron capture therapy (GdNCT) technique through a hybrid Monte Carlo (MC) simulation approach. At first, Snyder phantom including a spherical tumor was simulated by Geant4 MC code and relevant energy electron spectra to different 157Gd concentrations including 100, 250, 500, and 1000 ppm were calculated following the neutron irradiation of simulated phantom. Scored energy electron spectra were then imported to Monte Carlo damage simulation (MCDS) code to estimate RBE values (both RBESSB and RBEDSB) at different gadolinium concentrations and oxygen levels from 10 to 100%. The results indicate that variations of 157Gd can affect the energy spectrum of released secondary electrons including Auger electrons. Variation of gadolinium concentration from 100 to 1000 ppm in tumor region can change RBESSB and RBEDSB values by about 0.1% and 0.5%, respectively. Besides, maximum variations of 4.3% and 2% were calculated for RBEDSB and RBESSB when cell oxygen level changed from 10 to 100%. From the results, variations of considered gadolinium and oxygen concentrations during GdNCT can influence RBE values. Nevertheless, due to the not remarkable changes in the intensity of Auger electrons, a slight difference in RBE values would be expected at various 157Gd concentrations, although considerable RBE changes were calculated relevant to the oxygen alternations inside tumor tissue.

The dirty and clean dose concept: Towards creating proton therapy treatment plans with a photon-like dose response.

BACKGROUND: Applying tolerance doses for organs at risk (OAR) from photon therapy introduces uncertainties in proton therapy when assuming a constant relative biological effectiveness (RBE) of 1.1. PURPOSE: This work introduces the novel dirty and clean dose concept, which allows for creating treatment plans with a more photon-like dose response for OAR and, thus, less uncertainties when applying photon-based tolerance doses. METHODS: The concept divides the 1.1-weighted dose distribution into two parts: the clean and the dirty dose. The clean and dirty dose are deposited by protons with a linear energy transfer (LET) below and above a set LET threshold, respectively. For the former, a photon-like dose response is assumed, while for the latter, the RBE might exceed 1.1. To reduce the dirty dose in OAR, a MaxDirtyDose objective was added in treatment plan optimization. It requires setting two parameters: LET threshold and max dirty dose level. A simple geometry consisting of one target volume and one OAR in water was used to study the reduction in dirty dose in the OAR depending on the choice of the two MaxDirtyDose objective parameters during plan optimization. The best performing parameter combinations were used to create multiple dirty dose optimized (DDopt) treatment plans for two cranial patient cases. For each DDopt plan, 1.1-weighted dose, variable RBE-weighted dose using the Wedenberg RBE model and dose-average LETd distributions as well as resulting normal tissue complication probability (NTCP) values were calculated and compared to the reference plan (RefPlan) without MaxDirtyDose objectives. RESULTS: In the water phantom studies, LET thresholds between 1.5 and 2.5 keV/µm yielded the best plans and were subsequently used. For the patient cases, nearly all DDopt plans led to a reduced Wedenberg dose in critical OAR. This reduction resulted from an LET reduction and translated into an NTCP reduction of up to 19 percentage points compared to the RefPlan. The 1.1-weighted dose in the OARs was slightly increased (patient 1: 0.45 Gy(RBE), patient 2: 0.08 Gy(RBE)), but never exceeded clinical tolerance doses. Additionally, slightly increased 1.1-weighted dose in healthy brain tissue was observed (patient 1: 0.81 Gy(RBE), patient 2: 0.53 Gy(RBE)). The variation of NTCP values due to variation of α/ß from 2 to 3 Gy was much smaller for DDopt (2 percentage points (pp)) than for RefPlans (5 pp). CONCLUSIONS: The novel dirty and clean dose concept allows for creating biologically more robust proton treatment plans with a more photon-like dose response. The reduced uncertainties in RBE can, therefore, mitigate uncertainties introduced by using photon-based tolerance doses for OAR.